Resonance Offset Tailored Pulses for NMR Quantum Computation 



H. K. Cummins 

Oxford Centre for Quantum Computation, Clarendon Laboratory, Parks Road, Oxford 0X1 3PU, UK 

E-mail: h.cummms@physics.ox.ac.uk 



J. A. Jones^ 

Oxford Centre for Quantum Computation, Clarendon Laboratory, Parks Road, Oxford 0X1 3PU, UK, and 
Oxford Centre for Molecular Sciences, New Chemistry Laboratory, South Parks Road. Oxford, 0X1 3QT, UK 

E-mail: jonathan.joncs@qubit.org 



We describe novel composite pulse sequences which 
act as general rotors and thus are suitable for nuclear 
magnetic resonance (NMR) quantum computation. The 
Resonance Offset Tailoring To Enhance Nutations ap- 
proach permits perfect compensation of off-resonance 
errors at two selected frequencies placed symmetrically 
around the frequency of the RF source. 
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1. INTRODUCTION 

Composite pulses [1, 2] play an important role in many 
NMR experiments, as they allow the effects of experi- 
mental imperfections, such as pulse length errors and off- 
resonance effects, to be reduced. Such pulses can also 
prove useful in NMR implementations of quantum infor- 
mation processing devices, such as simple NMR quantum 
computers [3, 4, 5, 6, 7], where they act to reduce system- 
atic errors in quantum logic gates [8] . Unfortunately many 
conventional composite pulse sequences are not appropri- 
ate for quantum computers as they only perform well for 
certain initial states, while pulse sequences designed for 
quantum information processing must act as general ro- 
tors, that is they must perform well for any initial state. 

Composite pulses of this kind, which are sometimes 
called Class A composite pulses [1], are rarely needed for 
conventional NMR experiments, and so relatively little is 
known about them. One important example is a com- 
posite 90° pulse developed by Tycko [1, 9], which has re- 
cently been generalised to arbitrary rotation angles [8]. 
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These composite pulses give excellent compensation of off- 
resonance effects at small offset frequencies, such as those 
found for nuclei, but arc of no use for the much larger 
off resonance frequencies typically found for ^^C. 

Fortunately when composite pulses are used for NMR 
quantum computation one great simplification can be 
made: it is only necessary that the pulse sequence perform 
well over a small number of discrete frequency ranges, cor- 
responding to the resonance frequencies of the nuclei used 
to implement qubits; it is not necessary to design pulses 
which work well over the whole frequency range. In partic- 
ular it is quite common in NMR quantum computation to 
use at most two spins of each nuclear species (see, for ex- 
ample, [10]), and it is convenient to place the RF frequency 
in the centre of the spectrum, so that the two spins have 
equal and opposite resonance offsets [11]. Thus it suffices 
to tailor the composite pulse sequence to work well at these 
two frequencies; the performance at all other frequencies 
can be completely ignored. 

Here we explain how Resonance Offset Tailoring To 
Enhance Nutations may be used to produce composite 
pulse sequences which give perfect compensation of off- 
resonance effects. These ROTTEN pulses act as perfect 
general rotors at two frequencies, offset from the RF fre- 
quency by ±(5, and are well suited to NMR quantum com- 
putation; in combination with periods of free precession 
they provide an adequate set of gates, permitting any op- 
eration to be performed. ROTTEN pulses are simple to 
implement, and may be derived for any desired resonance 
offset as long as (5 < ^/2>Vl. 

2. RESULTS 

We choose to implement our composite rotation using 
a sequence of three radio-frequency (RF) pulses. In the 
absence of off-resonance effects any such pulse sequence 



1 



2 



CUMMINS AND JONES 



can be described by stating two angles describing each 
pulse, 9j, the nutation angle of the jth pulse, and (pj, the 
phase angle of the nutation axis in the xy-p\a,ne. In the 
presence of off-resonance effects it is convenient to retain 
this description, except that Oj is now a nominal nutation 
angle, and the nutation axis is no longer in the xy-plane 
(although remains a good description of the phase an- 
gle within the plane). It is also necessary to characterise 
the off-resonance behaviour, which is conveniently param- 
eterised using either the off-resonance fraction / = S/ui 
(where 5 is the off-resonance frequency, and i>i the nuta- 
tion rate, both measured in Hertz) or, equivalently, the RF 
tilt angle (the tilt angle of the nutation axis away from the 
2:-axis), given by tan(A) = I'l/S = 1/f. The propagator 
describing a single RF pulse is then 

(where Ix , ly and Iz and the conventional one spin product 
operators [12]), while the overall propagator describing the 
three pulse sequence is U = U3U2U1. 

In order to produce a perfectly compensated pulse it is 
necessary to find values of the six angles (^1, (pi, O2, <^2, 
03, and (j)^) such that U implements the desired rotation 
for the desired value of /. A general search over these six 
values would be a major task, but fortunately the prob- 
lem can be substantially simplified. The requirement that 
the composite pulse has identical effects on resonances at 
frequencies ±S, so that U{f) = U{—f), imposes major re- 
strictions on the allowed values. Furthermore, a family of 
solutions exists for which the first and last pulses are iden- 
tical, thus reducing the underlying search space to four 
independent values. Examining the form of the asymmet- 
ric response term U{f) — U{—f) suggests the choices 

ei = tt/vT+T^ (2) 

and 

cos(<^i - 02) = (1 - (3) 

These values can then be inserted back into the expression 
for U, the result equated with the desired propagator (ne- 
glecting any irrelevant overall phase term), and the equa- 
tions then solved for 62 and It is simplest to begin by 
solving for a 6x pulse, that is an ideal pulse with nutation 
angle 9 and phase angle 0; pulses with any other phase 
angle can then be created by simply adding a phase shift 
to all three pulses. This gives the results 

62 = e/^/iTp (4) 

and 



= ibarccos ^ J . (5) 

(Throughout this paper we will use the positive solution 
of equation 5). Finally combining equations 3 and 5 gives 

(j)2 = n - (6) 



Examining equations 3 and 5 reveals a limitation to this 
approach: the phase angles only have real solutions when 
I/I < a/S- It is, of course, possible that other compos- 
ite pulse families exist with larger ranges of applicability, 
but clearly there is a limit beyond which any three pulse 
sequence will cease to function. However, even the lim- 
ited range of / values described here is far greater than 
anything which can be achieved with conventional (non- 
tailored) composite pulses, and is likely to be adequate for 
most purposes. 

It must, however, be remembered that ROTTEN com- 
posite pulses are only effective at the resonance offsets 
for which they have been tailored; at other resonance off- 
sets these pulses may perform very poorly indeed. This is 
shown in figure 1 which plots the fidelity of a ROTTEN 
composite pulse sequence optimised for / = ±-\/3 over a 
range of values of /. As ROTTEN pulses are designed 
to act as general rotors it is necessary to use a fidelity 
measure which applies over all possible starting values; we 
have chosen to use the rotor fidelity, A, defined by Levitt 
[1]. The rotor fidelity of the simple pulse is low except at 
small values of /, while that for the composite pulse is high 
in two regions around ±\/3. This is ideal for some imple- 
mentations of NMR quantum computation, but unlikely 
to prove useful in most conventional NMR experiments. 
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FIG. 1. The fidelity, A, of a simple 90° pulse (dashed line) and 
a ROTTEN composite 90° pulse for a range of values of the off- 
resonance fraction, /; the ROTTEN pulse was tailored for the values 
/ = ±x/3. 



For the remainder of this paper we will consider compos- 
ite pulses tailored for the case f = ±V3; this is not only 
the limit of our approach (and so the case where ROT- 
TEN pulses give the greatest improvement in comparison 
with conventional pulses), but also a choice which results 
in a particularly simple sequence. To achieve an ideal 9^ 
rotation the values required are 9i — 037r/2, 62 — 9/2, 
4'i = 4''i4> ^-nd (1)2 — (j) + n. The operation of simple and 
ROTTEN composite 90° pulses arc shown by magnetiza- 
tion trajectories in figure 2 for initial states of Ix, ly and 
Iz- The magnetization trajectories are complicated, but 
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FIG. 2. Grapefruit plots showing magnetization trajectories for 90j pulses using simple pulses and ROTTEN composite pulses optimised 
for an off-resonance fraction / = simple pulses (a, b, c) with no ofF-resonance effects (/ = 0); simple pulses (d, e, f) in the presence 
of large positive off-resonance effects (/ = V3); ROTTEN pulses (g, h, i) in the presence of large positive off-resonance effects (/ = x/3); 
ROTTEN pulses (j, k, 1) in the presence of large negative off-resonance effects (/ = — -v/S). Initial states are Ix (a, d, g, j), (b, e, h, k) 
and Iz (c, f, i, 1). Note that for a perfect 90J pulse applied to Ix (a) the magnetization does not leave the a;-axis. 
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FIG. 3. Experimental ^^C spectra of ^^C labeled glycine in a home built 500 MHz (^H frequency) NMR spectrometer at the Oxford 
Centre for Molecular Sciences: (a) using a simple 90° pulse; (b) using a ROTTEN pulse. The frequency separation between the Ca and C 
multiplets was 18480 Hz and the RF pulse power was reduced to about 23% of its maximum value (around 13.5 kHz) so that the off-resonance 
fraction was / v^- 



those for ROTTEN composite pulses terminate in the cor- 
rect locations, while simple 90° pulses give extremely poor 
results at such large resonance offsets. 

Finally we show the performance of these pulse se- 
quences in an actual NMR experiment. Figure 3 shows 
^^C spectra of ^"^C labeled glycine acquired using simple 
and ROTTEN composite 90° excitation pulses. In order to 
emphasize the performance of our composite pulses these 
spectra were acquired with the RF pulse power reduced 
so that the off-resonance fraction was v^. Under these 
circumstances a simple 90° excitation results in phase er- 
rors of ±90°, while the ROTTEN pulse should give perfect 
compensation. In practice small phase errors are observed 
in the ROTTEN spectra; calculations suggest that these 
arise from pulse length errors. 



3. CONCLUSIONS 

Resonance offset tailoring provides a simple and effective 
approach for removing large off-resonance effects in sys- 
tems where NMR resonances occur at two well separated 
frequencies. Unlike most conventional composite pulses, 
ROTTEN pulses provide theoretically perfect compensa- 
tion, and act as perfect rotors. These properties make 
them well suited to NMR quantum computation. They 
are, however, unlikely to be useful in more conventional 
NMR experiments. 
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